Novel Methodology for Measuring Regional Myocardial Efficiency.

Our approach differs from the usual global measure of cardiac efficiency by using PET/MRI to measure efficiency of small pieces of cardiac tissue whose limiting size is equal to the spatial resolution of the PET scanner. We initiated a dynamic cardiac PET study immediately prior to the injection of 15.1 mCi of 11C-acetate acquiring data for 25 minutes while simultaneously acquiring MRI cine data. 1) A 3D finite element (FE) biomechanical model of the imaged heart was constructed by utilizing nonrigid deformable image registration to alter the Dassault Systèmes FE Living Heart Model (LHM) to fit the geometry in the cardiac MRI cine data. The patient specific FE cardiac model with estimates of stress, strain, and work was transformed into PET/MRI format. 2) A 1-tissue compartment model was used to calculate wash-in (K1) and the linear portion of the decay in the PET 11C-acetate time activity curve (TAC) was used to calculate the wash-out k2(mono) rate constant. K1 was used to calculate blood flow and k2(mono) was used to calculate myocardial volume oxygen consumption ( MVO2 ). 3) Estimates of stress and strain were used to calculate Myocardial Equivalent Minute Work ( MEMW ) and Cardiac Efficiency = MEMW/MVO2 was then calculated for 17 tissue segments of the left ventricle. The global MBF was 0.96 ± 0.15 ml/min/gm and MVO2 ranged from 8 to 17 ml/100gm/min. Six central slices of the MRI cine data provided a range of MEMW of 0.1 to 0.4 joules/gm/min and a range of Cardiac Efficiency of 6 to 18%.

Multiresolution spatiotemporal mechanical model of the heart as a prior to constrain the solution for 4D models of the heart.

In several nuclear cardiac imaging applications (SPECT and PET), images are formed by reconstructing tomographic data using an iterative reconstruction algorithm with corrections for physical factors involved in the imaging detection process and with corrections for cardiac and respiratory motion. The physical factors are modeled as coefficients in the matrix of a system of linear equations and include attenuation, scatter, and spatially varying geometric response. The solution to the tomographic problem involves solving the inverse of this system matrix. This requires the design of an iterative reconstruction algorithm with a statistical model that best fits the data acquisition. The most appropriate model is based on a Poisson distribution. Using Bayes Theorem, an iterative reconstruction algorithm is designed to determine the maximum a posteriori estimate of the reconstructed image with constraints that maximizes the Bayesian likelihood function for the Poisson statistical model. The a priori distribution is formulated as the joint entropy (JE) to measure the similarity between the gated cardiac PET image and the cardiac MRI cine image modeled as a FE mechanical model. The developed algorithm shows the potential of using a FE mechanical model of the heart derived from a cardiac MRI cine scan to constrain solutions of gated cardiac PET images.

Deformation Microscopy for Dynamic Intracellular and Intranuclear Mapping of Mechanics with High Spatiotemporal Resolution.

Structural heterogeneity is a hallmark of living cells that drives local mechanical properties and dynamic cellular responses. However, the robust quantification of intracellular mechanics is lacking from conventional methods. Here, we describe the development of deformation microscopy, which leverages conventional imaging and an automated hyperelastic warping algorithm to investigate strain history, deformation dynamics, and changes in structural heterogeneity within the interior of cells and cell nuclei. Using deformation microscopy, we found that partial or complete disruption of LINC complexes in cardiomyocytes in vitro and lamin A/C deficiency in myocytes in vivo abrogate dominant tensile loading in the nuclear interior. We also found that cells cultured on stiff substrates or in hyperosmotic conditions displayed abnormal strain burden and asymmetries at interchromatin regions, which are associated with active transcription. Deformation microscopy represents a foundational approach toward intracellular elastography, with the potential utility to provide mechanistic and quantitative insights in diverse mechanobiological applications.

Incorporation of the Living Heart Model into the 4D XCAT Phantom for Cardiac Imaging Research.

The 4D extended cardiac-torso (XCAT) phantom has provided a valuable tool to study the effects of anatomy and motion on medical images, especially cardiac motion. One limitation of the XCAT was that it did not have a physiological basis which to realistically simulate variations in cardiac function. In this work, we incorporate into the XCAT anatomy the four-chamber FE Living Heart Model (LHM) developed by the Living Heart Project (LHP). The LHM represents the state of the art in cardiac FE simulation because of its ability to accurately replicate the biomechanical motion of the entire heart and its variations. We create a new series of 4D phantoms capable of simulating patients with varying body sizes and shapes; cardiac positions, orientations, and dynamics. While extendable to other imaging modalities and technologies, our goal is to use the FE-enhanced XCAT models to investigate the optimal use of computed tomography (CT) for the evaluation of coronary artery disease (CAD). With the ability to simulate realistic, predictive, patient quality 4D imaging data, the enhanced XCAT models will enable optimization studies to identify the most promising systems or system parameters for further clinical validation.

In Vivo Multiscale and Spatially-Dependent Biomechanics Reveals Differential Strain Transfer Hierarchy in Skeletal Muscle.

Biological tissues have a complex hierarchical architecture that spans organ to subcellular scales and comprises interconnected biophysical and biochemical machinery. Mechanotransduction, gene regulation, gene protection, and structure-function relationships in tissues depend on how force and strain are modulated from macro to micro scales, and vice versa. Traditionally, computational and experimental techniques have been used in common model systems (e.g., embryos) and simple strain measures were applied. But the hierarchical transfer of mechanical parameters like strain in mammalian systems is largely unexplored in vivo. Here, we experimentally probed complex strain transfer processes in mammalian skeletal muscle tissue over multiple biological scales using complementary in vivo ultrasound and optical imaging approaches. An iterative hyperelastic warping technique quantified the spatially-dependent strain distributions in tissue, matrix, and subcellular (nuclear) structures, and revealed a surprising increase in strain magnitude and heterogeneity in active muscle as the spatial scale also increased. The multiscale strain heterogeneity indicates tight regulation of mechanical signals to the nuclei of individual cells in active muscle, and an emergent behavior appearing at larger (e.g. tissue) scales characterized by dramatically increased strain complexity.

The direct incorporation of perfusion defect information to define ischemia and infarction in a finite element model of the left ventricle.

This paper describes the process in which complex lesion geometries (specified by computer generated perfusion defects) are incorporated in the description of nonlinear finite element (FE) mechanical models used for specifying the motion of the left ventricle (LV) in the 4D extended cardiac torso (XCAT) phantom to simulate gated cardiac image data. An image interrogation process was developed to define the elements in the LV mesh as ischemic or infarcted based upon the values of sampled intensity levels of the perfusion maps. The intensity values were determined for each of the interior integration points of every element of the FE mesh. The average element intensity levels were then determined. The elements with average intensity values below a user-controlled threshold were defined as ischemic or infarcted depending upon the model being defined. For the infarction model cases, the thresholding and interrogation process were repeated in order to define a border zone (BZ) surrounding the infarction. This methodology was evaluated using perfusion maps created by the perfusion cardiac-torso (PCAT) phantom an extension of the 4D XCAT phantom. The PCAT was used to create 3D perfusion maps representing 90% occlusions at four locations (left anterior descending (LAD) segments 6 and 9, left circumflex (LCX) segment 11, right coronary artery (RCA) segment 1) in the coronary tree. The volumes and shapes of the defects defined in the FE mechanical models were compared with perfusion maps produced by the PCAT. The models were incorporated into the XCAT phantom. The ischemia models had reduced stroke volume (SV) by 18-59 ml. and ejection fraction (EF) values by 14-50% points compared to the normal models. The infarction models, had less reductions in SV and EF, 17-54 ml. and 14-45% points, respectively. The volumes of the ischemic/infarcted regions of the models were nearly identical to those volumes obtained from the perfusion images and were highly correlated (R² = 0.99).

Direct measurement of intranuclear strain distributions and RNA synthesis in single cells embedded within native tissue.

Nuclear structure and mechanics play a critical role in diverse cellular functions, such as organizing direct access of chromatin to transcriptional regulators. Here, we use a new, to our knowledge, hybrid method, based on microscopy and hyperelastic warping, to determine three-dimensional strain distributions inside the nuclei of single living cells embedded within their native extracellular matrix. During physiologically relevant mechanical loading to tissue samples, strain was transferred to individual nuclei, resulting in submicron distributions of displacements, with compressive and tensile strain patterns approaching a fivefold magnitude increase in some locations compared to tissue-scale stimuli. Moreover, nascent RNA synthesis was observed in the interchromatin regions of the cells studied and spatially corresponded to strain patterns. Our ability to measure large strains in the interchromatin space, which reveals that movement of chromatin in the nucleus may not be due to random or biochemical mechanisms alone, but may result from the transfer of mechanical force applied at a distant tissue surface.

A Comparison of Hyperelastic Warping of PET Images with Tagged MRI for the Analysis of Cardiac Deformation.

The objectives of the following research were to evaluate the utility of a deformable image registration technique known as hyperelastic warping for the measurement of local strains in the left ventricle through the analysis of clinical, gated PET image datasets. Two normal human male subjects were sequentially imaged with PET and tagged MRI imaging. Strain predictions were made for systolic contraction using warping analyses of the PET images and HARP based strain analyses of the MRI images. Coefficient of determination R (2) values were computed for the comparison of circumferential and radial strain predictions produced by each methodology. There was good correspondence between the methodologies, with R (2) values of 0.78 for the radial strains of both hearts and from an R (2) = 0.81 and R (2) = 0.83 for the circumferential strains. The strain predictions were not statistically different (P ≤ 0.01). A series of sensitivity results indicated that the methodology was relatively insensitive to alterations in image intensity, random image noise, and alterations in fiber structure. This study demonstrated that warping was able to provide strain predictions of systolic contraction of the LV consistent with those provided by tagged MRI Warping.

Generation of a suite of 3D computer-generated breast phantoms from a limited set of human subject data.

PURPOSE: The authors previously reported on a three-dimensional computer-generated breast phantom, based on empirical human image data, including a realistic finite-element based compression model that was capable of simulating multimodality imaging data. The computerized breast phantoms are a hybrid of two phantom generation techniques, combining empirical breast CT (bCT) data with flexible computer graphics techniques. However, to date, these phantoms have been based on single human subjects. In this paper, the authors report on a new method to generate multiple phantoms, simulating additional subjects from the limited set of original dedicated breast CT data. The authors developed an image morphing technique to construct new phantoms by gradually transitioning between two human subject datasets, with the potential to generate hundreds of additional pseudoindependent phantoms from the limited bCT cases. The authors conducted a preliminary subjective assessment with a limited number of observers (n = 4) to illustrate how realistic the simulated images generated with the pseudoindependent phantoms appeared. METHODS: Several mesh-based geometric transformations were developed to generate distorted breast datasets from the original human subject data. Segmented bCT data from two different human subjects were used as the "base" and "target" for morphing. Several combinations of transformations were applied to morph between the "base' and "target" datasets such as changing the breast shape, rotating the glandular data, and changing the distribution of the glandular tissue. Following the morphing, regions of skin and fat were assigned to the morphed dataset in order to appropriately assign mechanical properties during the compression simulation. The resulting morphed breast was compressed using a finite element algorithm and simulated mammograms were generated using techniques described previously. Sixty-two simulated mammograms, generated from morphing three human subject datasets, were used in a preliminary observer evaluation where four board certified breast radiologists with varying amounts of experience ranked the level of realism (from 1 = "fake" to 10 = "real") of the simulated images. RESULTS: The morphing technique was able to successfully generate new and unique morphed datasets from the original human subject data. The radiologists evaluated the realism of simulated mammograms generated from the morphed and unmorphed human subject datasets and scored the realism with an average ranking of 5.87 ± 1.99, confirming that overall the phantom image datasets appeared more "real" than "fake." Moreover, there was not a significant difference (p > 0.1) between the realism of the unmorphed datasets (6.0 ± 1.95) compared to the morphed datasets (5.86 ± 1.99). Three of the four observers had overall average rankings of 6.89 ± 0.89, 6.9 ± 1.24, 6.76 ± 1.22, whereas the fourth observer ranked them noticeably lower at 2.94 ± 0.7. CONCLUSIONS: This work presents a technique that can be used to generate a suite of realistic computerized breast phantoms from a limited number of human subjects. This suite of flexible breast phantoms can be used for multimodality imaging research to provide a known truth while concurrently producing realistic simulated imaging data.

An analysis of the mechanical parameters used for finite element compression of a high-resolution 3D breast phantom.

PURPOSE: The authors previously introduced a methodology to generate a realistic three-dimensional (3D), high-resolution, computer-simulated breast phantom based on empirical data. One of the key components of such a phantom is that it provides a means to produce a realistic simulation of clinical breast compression. In the current study, they have evaluated a finite element (FE) model of compression and have demonstrated the effect of a variety of mechanical properties on the model using a dense mesh generated from empirical breast data. While several groups have demonstrated an effective compression simulation with lower density finite element meshes, the presented study offers a mesh density that is able to model the morphology of the inner breast structures more realistically than lower density meshes. This approach may prove beneficial for multimodality breast imaging research, since it provides a high level of anatomical detail throughout the simulation study. METHODS: In this paper, the authors describe methods to improve the high-resolution performance of a FE compression model. In order to create the compressible breast phantom, dedicated breast CT data was segmented and a mesh was generated with 4-noded tetrahedral elements. Using an explicit FE solver to simulate breast compression, several properties were analyzed to evaluate their effect on the compression model including: mesh density, element type, density, and stiffness of various tissue types, friction between the skin and the compression plates, and breast density. Following compression, a simulated projection was generated to demonstrate the ability of the compressible breast phantom to produce realistic simulated mammographic images. RESULTS: Small alterations in the properties of the breast model can change the final distribution of the tissue under compression by more than 1 cm; which ultimately results in different representations of the breast model in the simulated images. The model properties that impact displacement the most are mesh density, friction between the skin and the plates, and the relative stiffness of the different tissue types. CONCLUSIONS: The authors have developed a 3D, FE breast model that can yield high spatial resolution breast deformations under uniaxial compression for imaging research purposes and demonstrated that small changes in the mechanical properties can affect images generated using the phantom.

Incorporation of a left ventricle finite element model defining infarction into the XCAT imaging phantom.

The 4D extended cardiac-torso (XCAT) phantom was developed to provide a realistic and flexible model of the human anatomy and cardiac and respiratory motions for use in medical imaging research. A prior limitation to the phantom was that it did not accurately simulate altered functions of the heart that result from cardiac pathologies such as coronary artery disease (CAD). We overcame this limitation in a previous study by combining the phantom with a finite-element (FE) mechanical model of the left ventricle (LV) capable of more realistically simulating regional defects caused by ischemia. In the present work, we extend this model giving it the ability to accurately simulate motion abnormalities caused by myocardial infarction (MI), a far more complex situation in terms of altered mechanics compared with the modeling of acute ischemia. The FE model geometry is based on high resolution CT images of a normal male subject. An anterior region was defined as infarcted and the material properties and fiber distribution were altered, according to the bio-physiological properties of two types of infarction, i.e., fibrous and remodeled infarction (30% thinner wall than fibrous case). Compared with the original, surface-based 4D beating heart model of the XCAT, where regional abnormalities are modeled by simply scaling down the motion in those regions, the FE model was found to provide a more accurate representation of the abnormal motion of the LV due to the effects of fibrous infarction as well as depicting the motion of remodeled infarction. In particular, the FE models allow for the accurate depiction of dyskinetic motion. The average circumferential strain results were found to be consistent with measured dyskinetic experimental results. Combined with the 4D XCAT phantom, the FE model can be used to produce realistic multimodality sets of imaging data from a variety of patients in which the normal or abnormal cardiac function is accurately represented.

Realistic Simulation of Regional Myocardial Perfusion Defects for Cardiac SPECT Studies.

The current 3D XCAT phantom allows users to manually define the regional myocardial perfusion defect (MPD) as a simple pie-shaped wedge region with reduced activity level in the myocardium of left ventricle. To more accurately and realistically model the MPD, we have developed a new regional MPD model for the 3D XCAT phantom for myocardial perfusion SPECT (MP-SPECT) studies based on the location and the severity of the stenosis in a computer generated coronary arterial tree. First, we generated a detailed coronary arterial tree by extending the large proximal branches segmented from patient CT images to cover the whole heart using an iterative rule-based algorithm. Second, we determined the affected downstream vascular segments of a given stenosis. Third, we computed the activity of each myocardial region as a function of the inverse-distance-weighted average of the flow of the neighboring vascular segments. Fourth, we generated a series of bull's-eye maps of MP-SPECT images of different coronary artery stenosis scenarios. Fifth, we had expert physician readers to qualitatively assess the bull's-eye maps based on their similarity to typical clinical cases in terms of the shape, the extent, and the severity of the MPDs. Their input was used to iteratively revise the coronary artery tree model so that the MPDs were closely matched to those found in bull's-eye maps from patient studies. Finally, from our simulated MP-SPECT images, we observed that (1) the locations of the MPDs caused by stenoses at different main arteries were different largely according to their vascular territories, (2) a stenosis at a proximal branch produced a larger MPD than the one at a distal branch, and (3) a more severe stenosis produced a larger MPD than the less severe one. These observations were consistent to those found in clinical cases. Therefore, this new regional MPD model has enhanced the generation of realistic pathological MP-SPECT images using the XCAT phantom. When combining with the mechanical model of the myocardium, the new model can be extended for the simulation of 4D gated MP-SPECT simulation of a pathological heart with both perfusion and motion defects.

Strain measurement in the left ventricle during systole with deformable image registration.

The objective of this study was to validate a deformable image registration technique, termed Hyperelastic Warping, for left ventricular strain measurement during systole using cine-gated, non-tagged MR images with strains measured from tagged MRI. The technique combines deformation from high resolution, non-tagged MR image data with a detailed computational model, including estimated myocardial material properties, fiber direction, and active fiber contraction, to provide a comprehensive description of myocardial contractile function. A normal volunteer (male, age 30) with no history of cardiac pathology was imaged with a 1.5 T Siemens Avanto clinical scanner using a TrueFISP imaging sequence and a 32-channel cardiac coil. Both tagged and non-tagged cine MR images were obtained. The Hyperelastic Warping solution was evolved using a series of non-tagged images in ten intermediate phases from end-diastole to end-systole. The solution may be considered as ten separate warping problems with multiple templates and targets. At each stage, an active contraction was initially applied to a finite element model, and then image-based warping penalty forces were utilized to generate the final registration. Warping results for circumferential strain (R(2)=0.75) and radial strain (R(2)=0.78) were strongly correlated with results obtained from tagged MR images analyzed with a Harmonic Phase (HARP) algorithm. Results for fiber stretch, LV twist, and transmural strain distributions were in good agreement with experimental values in the literature. In conclusion, Hyperelastic Warping provides a unique alternative for quantifying regional LV deformation during systole without the need for tags.

Measuring regional changes in the diastolic deformation of the left ventricle of SHR rats using microPET technology and hyperelastic warping.

The objective of this research was to assess applicability of a technique known as hyperelastic warping for the measurement of local strains in the left ventricle (LV) directly from microPET image data sets. The technique uses differences in image intensities between template (reference) and target (loaded) image data sets to generate a body force that deforms a finite element (FE) representation of the template so that it registers with the target images. For validation, the template image was defined as the end-systolic microPET image data set from a Wistar Kyoto (WKY) rat. The target image was created by mapping the template image using the deformation results obtained from a FE model of diastolic filling. Regression analysis revealed highly significant correlations between the simulated forward FE solution and image derived warping predictions for fiber stretch (R (2) = 0.96), circumferential strain (R (2) = 0.96), radial strain (R (2) = 0.93), and longitudinal strain (R (2) = 0.76) (p < 0.001 for all cases). The technology was applied to microPET image data of two spontaneously hypertensive rats (SHR) and a WKY control. Regional analysis revealed that, the lateral freewall in the SHR subjects showed the greatest deformation compared with the other wall segments. This work indicates that warping can accurately predict the strain distributions during diastole from the analysis of microPET data sets.

Noninvasive determination of ligament strain with deformable image registration.

Ligament function and propensity for injury are directly related to regional stresses and strains. However, noninvasive techniques for measurement of strain are currently limited. This study validated the use of Hyperelastic Warping, a deformable image registration technique, for noninvasive strain measurement in the human medial collateral ligament using direct comparisons with optical measurements. Hyperelastic Warping determines the deformation map that aligns consecutive images of a deforming material, allowing calculation of strain. Diffeomorphic deformations are ensured by representing the deformable image as a hyperelastic material. Ten cadaveric knees were subjected to six loading scenarios each. Tissue deformation was documented with magnetic resonance imaging (MRI) and video-based experimental measurements. MRI datasets were analyzed using Hyperelastic Warping, representing the medial collateral ligament (MCL) with a hexahedral finite element (FE) model projected to a manually segmented ligament surface. The material behavior was transversely isotropic hyperelastic. Warping predictions of fiber stretch were strongly correlated with experimentally measured strains (R (2) = 0.81). Both sets of measurements were in agreement with previous ex vivo studies. Warping predictions of fiber stretch were insensitive to bulk:shear modulus ratio, fiber stiffness, and shear modulus in the range of +2.5SD to -1.0SD. Correlations degraded when the shear modulus was decreased to 2.5SD below the mean (R (2) = 0.56), and when an isotropic constitutive model was substituted for the transversely isotropic model (R (2) = 0.65). MCL strains in the transitional region near the joint line, where the material behavior and material symmetry are more complex, showed the most sensitivity to changes in shear modulus. These results demonstrate that Hyperelastic Warping requires the use of a constitutive model that reflects the material symmetry, but not subject-specific material properties for accurate strain predictions for this application. Hyperelastic Warping represents a powerful technique for noninvasive strain measurement of musculoskeletal tissues and has many advantages over other image-based strain measurement techniques.

Normal and pathological NCAT image and phantom data based on physiologically realistic left ventricle finite-element models.

The four-dimensional (4-D) NURBS-based cardiac-torso (NCAT) phantom, which provides a realistic model of the normal human anatomy and cardiac and respiratory motions, is used in medical imaging research to evaluate and improve imaging devices and techniques, especially dynamic cardiac applications. One limitation of the phantom is that it lacks the ability to accurately simulate altered functions of the heart that result from cardiac pathologies such as coronary artery disease (CAD). The goal of this work was to enhance the 4-D NCAT phantom by incorporating a physiologically based, finite-element (FE) mechanical model of the left ventricle (LV) to simulate both normal and abnormal cardiac motions. The geometry of the FE mechanical model was based on gated high-resolution X-ray multislice computed tomography (MSCT) data of a healthy male subject. The myocardial wall was represented as a transversely isotropic hyperelastic material, with the fiber angle varying from -90 degrees at the epicardial surface, through 0 degrees at the midwall, to 90 degrees at the endocardial surface. A time-varying elastance model was used to simulate fiber contraction, and physiological intraventricular systolic pressure-time curves were applied to simulate the cardiac motion over the entire cardiac cycle. To demonstrate the ability of the FE mechanical model to accurately simulate the normal cardiac motion as well as the abnormal motions indicative of CAD, a normal case and two pathologic cases were simulated and analyzed. In the first pathologic model, a subendocardial anterior ischemic region was defined. A second model was created with a transmural ischemic region defined in the same location. The FE-based deformations were incorporated into the 4-D NCAT cardiac model through the control points that define the cardiac structures in the phantom which were set to move according to the predictions of the mechanical model. A simulation study was performed using the FE-NCAT combination to investigate how the differences in contractile function between the subendocardial and transmural infarcts manifest themselves in myocardial Single photon emission computed tomography (SPECT) images. The normal FE model produced strain distributions that were consistent with those reported in the literature and a motion consistent with that defined in the normal 4-D NCAT beating heart model based on tagged magnetic resonance imaging (MRI) data. The addition of a subendocardial ischemic region changed the average transmural circumferential strain from a contractile value of -0.09 to a tensile value of 0.02. The addition of a transmural ischemic region changed average circumferential strain to a value of 0.13, which is consistent with data reported in the literature. Model results demonstrated differences in contractile function between subendocardial and transmural infarcts and how these differences in function are documented in simulated myocardial SPECT images produced using the 4-D NCAT phantom. Compared with the original NCAT beating heart model, the FE mechanical model produced a more accurate simulation for the cardiac motion abnormalities. Such a model, when incorporated into the 4-D NCAT phantom, has great potential for use in cardiac imaging research. With its enhanced physiologically based cardiac model, the 4-D NCAT phantom can be used to simulate realistic, predictive imaging data of a patient population with varying whole-body anatomy and with varying healthy and diseased states of the heart that will provide a known truth from which to evaluate and improve existing and emerging 4-D imaging techniques used in the diagnosis of cardiac disease.

Measurement of strain in the left ventricle during diastole with cine-MRI and deformable image registration.

The assessment of regional heart wall motion (local strain) can localize ischemic myocardial disease, evaluate myocardial viability, and identify impaired cardiac function due to hypertrophic or dilated cardiomyopathies. The objectives of this research were to develop and validate a technique known as hyperelastic warping for the measurement of local strains in the left ventricle from clinical cine-magnetic resonance imaging (MRI) image datasets. The technique uses differences in image intensities between template (reference) and target (loaded) image datasets to generate a body force that deforms a finite element (FE) representation of the template so that it registers with the target image. To validate the technique, MRI image datasets representing two deformation states of a left ventricle were created such that the deformation map between the states represented in the images was known. A beginning diastolic cine-MRI image dataset from a normal human subject was defined as the template. A second image dataset (target) was created by mapping the template image using the deformation results obtained from a forward FE model of diastolic filling. Fiber stretch and strain predictions from hyperelastic warping showed good agreement with those of the forward solution (R2=0.67 stretch, R2=0.76 circumferential strain, R2=0.75 radial strain, and R2=0.70 in-plane shear). The technique had low sensitivity to changes in material parameters (deltaR2= -0.023 fiber stretch, deltaR2=-0.020 circumferential strain, deltaR2=-0.005 radial strain, and deltaR2=0.0125 shear strain with little or no change in rms error), with the exception of changes in bulk modulus of the material. The use of an isotropic hyperelastic constitutive model in the warping analyses degraded the predictions of fiber stretch. Results were unaffected by simulated noise down to a signal-to-noise ratio (SNR) of 4.0 (deltaR2= -0.032 fiber stretch, deltaR2=-0.023 circumferential strain, deltaR2=-0.04 radial strain, and deltaAR2=0.0211 shear strain with little or no increase in rms error). This study demonstrates that warping in conjunction with cine-MRI imaging can be used to determine local ventricular strains during diastole.

Strain measurement in coronary arteries using intravascular ultrasound and deformable images.

Atherosclerotic plaque rupture is responsible for the majority of myocardial infarctions and acute coronary syndromes. Rupture is initiated by mechanical failure of the plaque cap, and thus study of the deformation of the plaque in the artery can elucidate the events that lead to myocardial infarction. Intravascular ultrasound (IVUS) provides high resolution in vitro and in vivo cross-sectional images of blood vessels. To extract the deformation field from sequences of IVUS images, a registration process must be performed to correlate material points between image pairs. The objective of this study was to determine the efficacy of an image registration technique termed Warping to determine strains in plaques and coronary arteries from paired IVUS images representing two different states of deformation. The Warping technique uses pointwise differences in pixel intensities between image pairs to generate a distributed body force that acts to deform a finite element model. The strain distribution estimated by image-based Warping showed excellent agreement with a known forward finite element solution, representing the gold standard, from which the displaced image was created. The Warping technique had a low sensitivity to changes in material parameters or material model and had a low dependency on the noise present in the images. The Warping analysis was also able to produce accurate strain distributions when the constitutive model used for the Warping analysis and the forward analysis was different. The results of this study demonstrate that Warping in conjunction with in vivo IVUS imaging will determine the change in the strain distribution resulting from physiological loading and may be useful as a diagnostic tool for predicting the likelihood of plaque rupture through the determination of the relative stiffness of the plaque constituents.